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Understanding and designing inertial confinement fusion (ICF) implosions through radiation- 
hydrodynamics simulations rely on the accurate knowledge of the equation of state (EOS) of the 
deuterium and tritium fuels. To minimize the drive energy for ignition, the imploding shell of DT- 
fuel needs to be kept as cold as possible. Such low-adiabat ICF implosions can access to coupled and 
degenerate plasma conditions, in which the analytical or chemical EOS models become inaccurate. 
Using the path integral Monte Carlo (PIMC) simulations we have derived a first-principles EOS 
(FPEOS) table of deuterium that covers typical ICF fuel conditions at densities ranging from 0.002 
to 1596 g/cm'^ and temperatures of 1.35 eV — 5.5 keV. We report the internal energy and the 
pressure, and discuss the structure of the plasma in terms of pair correlation functions. When 
compared with the widely used SESAME table and the revised Kerley03 table, discrepancies in 
the internal energy and in the pressure are identified for moderately coupled and degenerate plasma 
conditions. In contrast to the SESAME table, the revised Kerley03 table is in better agreement with 
our FPEOS results over a wide range of densities and temperatures. Although subtle differences still 
exist for lower temperatures (T < 10 eV) and moderate densities (1 — 10 g/cm'^), hydrodynamics 
simulations of cryogenic ICF implosions using the FPEOS table and the KerleyOS table have resulted 
in similar results for the peak density, areal density pR, and neutron yield, which are significantly 
different from the SESAME simulations. 

PACS numbers: 52.25.Kn, 51.30.+i, 62.50.-p, 64.10. +h 



I. INTRODUCTION 

Inertial confinement fusion (ICF) has been pursued for 
decades since the concept was introduced in 1972^. In 
the traditional central-hot-spot ignition designs, a cap- 
sule of cryogenic deuterium-tritium (DT) covered with 
plastic ablator is driven to implode either directly by in- 
tense laser pulses^ or indirectly by x-rays in a hohlraum'^. 
To minimize the driving energy required for ignition, the 
imploding DT-capsule needs to be maintained as cold 
as possibl^ for high compressions (larger than a thou- 
sand times that of the solid DT density) at the stag- 
nation stage. This can either be done with fine-tuned 
shock^or with ramp compression waves. The reduction 
in temperature leads to pressures in the imploding DT- 
shell that are just above the Fermi degeneracy pressure. 
This is conventionally characterized by the so-called adia- 
bat parameter a = P/Pp- Low-adiabat ICF designs with 
1 < a < 2 are currently studied with indirect-drive im- 
plosions at the National Ignition Facility (NIFjP. Direct- 
drive ignition designJ^ for NIF also place the DT-shell 
adiabat at a low value of 2 < a < 3. Cryogenic DT 
targets scaled from the hydro-equivalent NIF designs are 
routinely imploded with a direct drive at the Omega laser 
facilitjM 

Since the compressibility of a material is determined 
by its equation of state (EOS)^, the accurate knowledge 
of the EOS of the DT-fuel is essential for designing ICF 
ignition targets and predicting the performance of the 
target during ICF implosions. To perform radiation- 
hydrodynamics simulations of ICF implosions, one needs 
to know the pressure and energy of the DT-fuel and the 



ablator materials at various density and temperature con- 
ditions, which are usually provided by EOS tables or an- 
alytical formulas. Various EOS tables for deuterium have 
been assembled because its importance in ICF applica- 
tions, planetary science and high pressure physics. 

The widely used SESAME EOS table of deuteriunPESl 
was based upon a chemical model for hydroger^^ that 
describes the material in terms of well-defined chemical 
species like II2 molecules, H atoms and free protons and 
electrons. Their interaction as well as many-body and 
degeneracy effects are treated approximately. For the 
SESAME table, liquid perturbation theory was adopted 
in the molecular/atomic fluid phase for ICF plasma con- 
ditions. A flrst-order expansion that only takes into ac- 
count nearest neighbor interactions was used in the orig- 
inal SESAME tablfP. 

Chemical models are expected to work well in the 
regime of weak coupling. However, in ICF implosions, 
the DT shell goes through a wide range of densities from 
0.1 up to 1000 g/cm^ and temperatures varying from a 
few electron volts (eV) to several hundreds of electron 
volt^^m^ which include plasma conditions with moder- 
ately strong coupling. This provides the primary mo- 
tivation for this paper, where we derive the deuterium 

from first-principles path integral Monte Carlo simula- 
tionpHH. 

The conditions for a low-adiabat {a ~ 2.5) cryogenic 
DT implosion on OMEGA are shown in Fig. [T] on pan- 
els (a)-(c). Panels (d)-(f) characterize the conditions 
for a direct-drive ignition design for NIF that is hydro- 
equivalent to the OMEGA implosion. In panels (a) and 
(d), we plot the laser pulse shapes. Panels (b) and (e) 
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FIG. 1: (Color online) (a)-(c): A cryogenic DT-implosion on 
OMEGA with the triple-picket step-pulse; (d)-(f): A direct- 
drive ignition design for NIF, scaled from hydro-equivalent 
OMEGA implosions. In both cases, strongly coupled and 
degenerated plasma conditions are indeed accessed. 



show the density, p, and temperature, T, path of the 
driven DT-sheh that we derived with one dimensional 
(ID) hydro-simulations using the hydro-code LILACP. 
The DT shell is predicted to undergo a variety of drive 
stages including several shocks and the final push by the 
main pulse. 

The p-T path of the imploding DT shell can be pro- 
jected onto a plane spanned by the coupling parameter 
r = l/{rskbT) and the degeneracy parameter 6 = T/Tp. 
Tp = X (STT^n)^/'^ is the Fermi temperature of the 

electrons in a fully ionized plasma and is the Wigner- 
Seitz radius that is related to the number density of the 
electrons, n — 3/{ATrr^). One finds that the imploding 
shells indeed pass through the strongly coupled (F > 1) 
and degenerate {6 < 1) regimes and one expects cou- 
pling and degeneracy effects to play significant roles in 
the compression and yield-production in low-adiabat ICF 
implosiond^. 

Strong coupling and degeneracy effects in ICF plasmas 
have recently attracted much attention, as they may re- 
define the so-called ID-physics of ICF implosions. The 
essential pieces of physics models used in ICF hydro- 
simulations, such as the electron-ion energy relaxation 
rat^^, the thermal conductivity'^, the fusion reaction 
ratJ^H^ and viscosity and mutual diffusion in deuterium- 
tritium mixtureJ^ in coupled and degenerated plasmas 
have been re-examined recently with experimental and 
theoretical methods. EOS measurements of liquid deu- 
terium along the principal Hugoniot reaching about 100- 
200 GPa have been performed using laser- driven shock 



wave^^imni^ m agneti cally driven flyer^^^^, and conver- 
gent explosivesS^ED, First-principles computer simula- 
tions have emerged as the preferred theoretical tool to 



derive the EOS of deuterium under such extreme con- 
ditions. Two methods have been most successful: den- 
sity functional molecular dynamics (DFT-MD^pSEll and 
the path integral Monte Carlo (PIMCj^^. In contrast 
to chemical models, these first-principles methods can 
take many-body effects fully into account. Results from 
such simulations have also been used to revise the original 
SESAME EOS table of deuterium to yield the improved 
Kerley03 EOS tabltP. 

For ICF applications, we are especially concerned 
about the EOS accuracy along the implosion path in 
the density-temperature plane, i.e., in the range of p = 
0. 1 - 1000 g/cm^ and T = 1 - 1000 eV. For such high tem- 
peratures, standard DFT methods become prohibitively 
expansive because of the large number of electronic or- 
bitals that would need to be included in the calcula- 
tion to account for electronic excitation^^. Orbital- 
free semi-classical simulation methods based on Thomas- 
Fermi theorjB^ is more efficient but they approximate 
electronic correlation effects and cannot represent chem- 
ical bonds. Therefore, in current form, they cannot de- 
scribe the systems at lower temperatures accurately. 

Path integral Monte Carlo has been shown to work 
rather well for EOS calcul ation s of low-Z materials such 
as deuteriunPS! and heliunP^12II_ jn this paper we present 
a first-principles equation of state (FPEOS) table of deu- 
terium from restricted PIMC simulations^. This method 
has been successfully applied to compute the deuterium 
EOgilMSI up to a density of p = 5.388 g/cm^. At lower 
temperatures, the PIMC results have been shown to 
agree well with DFT-MD calculations for hydrogerP^ and 
more recently for hclium^^. 

Our FPEOS table derived from PIMC covers the whole 
DT-shell plasma conditions throughout the low-adiabat 
ICF implosions. Specifically, our table covers densities 
ranging from 0.002 to 1596 gcm~^ and temperatures of 
1.35 eV — 5.5 keV. When compared with the widely used 
SESAME-EOS table and the revised Kerley03-EOS ta- 
ble, discrepancies in the internal energy and the pressure 
have been identified in moderately coupled and degen- 
erate regimes. Hydrodynamics simulations for cryogenic 
ICF implosions using our FPEOS table and the KerleyOS- 
EOS table have resulted in similar peak density, areal 
density pR, and neutron yield, which differ significantly 
from the SESAME simulations. 

The paper is organized as follows. A brief descrip- 
tion of the path integral Monte Carlo method is given in 
Sec. II. In Sec. Ill our FPEOS table is presented. In 
Sec. IV, we characterize the properties of the deuterium 
plasma for a variety of density and temperature condi- 
tions in terms of pair correlation functions. Comparisons 
between the FPEOS table, the SESAME and the Kcr- 
ley03 EOS as well as the simple Debye-Hiickel plasma 
model are made in Sec. V. In Sec. VI, we analyze the 
implications of different EOS tables for ICF applications 
through hydro-simulations and comparisons with exper- 
iments. The paper is summarized in Sec. VII. 
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II. THE PATH INTEGRAL MONTE CARLO 
METHOD 

Path integral Monte Carlo (PIMC) is the appropri- 
ate computational technique for simulating many-body 
quantum systems at finite temperatures. In PIMC calcu- 
lations, electrons and ions are treated on equal footing as 
paths, which means the quantum effects of both species 
are included consistently, although for the temperatures 
under consideration, the zero-point motion and exchange 
effects of the nuclei are negligible. 

The fundamental idea of the path integral approach is 
that the density matrix of a quantum system at temper- 
ature, T, can be expressed as a convolution of density 
matrices at a much higher temperature, M x T: 

p(R, R'; /?) = y" dtii dR.2 • • • dR-M-i p(R, Ri; A/J) 
X p(Ri,R2;A/3)...p(RM-i,R';A/3).(l) 

This is an exact expression. The integral on the right can 
be interpreted as a weighted average over all paths that 
connect the points R and R'. R is a collective variable 
that denote the positions of all particles R = {ri, . . . rjv}. 
j3 = l/khT represents length of the path in "imaginary 
time" and A/3 = /3/M is the size of each of the M time 
steps. 

From the free particle density matrix which can be 
used for the high-temperature density-matrices. 



exp 



./\2 



(2) 



one can estimate that the separation of two adjacent posi- 
tions on the path, Ar — ri+i — can only be on the order 

of J lS.l3/m while the separation of the two end points 



is approximately \Jh^l3/in. One can consequently inter- 
pret the positions Ri . . . Rj\,/_i as intermediate points on 
a path from R and R'. The multi-dimensional integra- 
tion over all paths in Eq. IT] c an be performed efhciently 
with Monte Carlo methoda^Zl. 

In general observables associated with operator, O, can 
be derived from. 



O 



J dR/ dR' (R| 6 |R') p(R', R; /3) 
/dR p(R, R; /3) 



(3) 



but for the kinetic and potential energies, Ex and Ep, as 
well as for pair correlation functions only diagonal matrix 
elements (R = R') are needed. The total internal energy 
follows from E ~ Ek + Ep and the pressure, P, can be 
obtained from the virial theorem for Coulomb systems. 



P = {2Ek + Ep)lW 



(4) 



Electrons are fermions and their fermionic characters 
matters for the degenerate plasma conditions under con- 
sideration. This implies one needs to construct an anti- 
symmetric many-body density matrix, which can be de- 
rived by introducing a sum of all permutations, V , and 
then also include paths from R to VTil . While this ap- 
proach works well for boson^^, for fermions each permu- 
tation must be weighted by a factor (—1)^. The partial 
cancellation of contributions with opposite signs leads to 
an extremely inefhcient algorithm when the combined po- 
sition and permutation space is sampled directly. This is 
known as Fermion sign problem, and its severity increases 
as the plasma becomes more degenerate. 

We deal with the Fermion sign p roblem by introducing 
the fixed node approximatiorPSHSl^ 

pp{R,R';f3) = ^ E (-1)'' / (5) 

PT(R,Rt;t)>0 

where one only includes those paths that satisfy the nodal 
constraint, /Ot(R, Ri;0 > 0, at every point. S'[Rt] is the 
action of the path and px is a fermionic trial density a 
matrix that must be given in analytic form. For this 
paper, we rely on free particle nodes. 



pt(R,R';/3) 



pW(ri,r'i;/3) ... (r^v, r'^; /3) 



W(ri,rW;/3) 



■ATI / 



(6) 



V is the volume. 



but the nodes of a variational density matrbP^l have also 
been employed in PIMC computationji^llllll. 

We have performed a number of convergence tests to 
minimize errors from using a finite time step and from 
a finite number of particles in cubic simulation cells 
with periodic boundary conditions. We determined a 
time step of A/3 < [100 x kbTp]"^ was sufficient to ac- 
curately account for all interactions and degeneracy ef- 
fects. We perform our PIMC calculations with different 
numbers of atoms depending on the deuterium density: 

= 64 atoms for p < 2.5 gcm~^, TV = 128 atoms for 
2.5 < p < 10.5 gcm"^. and N = 256 atoms for p > 10.5 
gcm^'^. 



III. THE FPEOS TABLE OF DEUTERIUM 

We have carried out PIMC calculations for a variety of 
density and temperature conditions that are of interests 
to inertial confinement fusion applications. The result- 
ing FPEOS table for deuterium covers the density range 
from 0.0019636 gcm^^^ (r^ = 14 in units of Bohr radii 
oo) to 1596.48802 gcm^^ ^ o.l5 Qq) and the tem- 
perature interval from 15 625 K (~ 1.35 eV) to 6.4 x 10'' 
K (~ 5515.09 eV). Fig. [2] shows the conditions for every 
simulation combined with lines for F = 1 and 9 = 1 to 
indicate the boundaries between coupled/uncoupled and 
degenerate/non-degenerate plasma conditions. Plasma 
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FIG. 2: (Color online) The temperature and density condi- 
tions covered by the FPEOS table. The gray circles repre- 
sent our PIMC calculations, while the shell conditions in ICF 
implosions are schematically shown by the region in orange 
color. The blue and green lines of 9 = 1 and F = 1 charac- 
terize the boundaries of degeneracy and coupling conditions, 
respectively. 



conditions in the upper left corner of the diagram are 
weakly coupled and classical (F ^ 1 and 9^1), while 
the lower right of the diagram represent strongly coupled 
and highly degenerate conditions (F ^ 1 and 6 ^ 1). 
The lowest temperatures in our PIMC calculations reach 
to the regime of ~ 0.1. 

To give an example for the ICF plasma conditions, we 
added the conditions of an imploding DT capsule shown 
in Fig. [l] to Fig. [2| It can be seen that the DT sheU 
undergoes a change from strongly coupled to an uncou- 
pled regime during the shock transits. The electronic 
conditions change from fully degenerated to partially de- 
generate accordingly. All these conditions are covered by 
our PIMC results that we have assembled into the follow- 
ing FPEOS table |lj The pressure and the internal energy 
as well as their statistical error bars from PIMC simu- 
lations are listed for different density and temperature 
conditions. 



IV. PARTICLE CORRELATIONS 

The correlation functions, g{r), between different pairs 
of particles such as electron-electron, electron-ion, and 
ion-ion are particularly interesting for analyzing the 
physical and chemical changes in the plasma at various 
density and temperature conditions. The g{r) are avail- 
able directly in PIMC simulations. We first show the 
density effects on the structure of the fluid structure by 
showing how the g{r) functions change with density for 
three temperatures of 15 625 K, 2.5 x 10^ K, and 2 x 10^ 
K in Figs. |3][5l 

Fig. |3ja) shows a clear peak in the ion-ion correlation 
function, gddii"), for 0.1 g/cnv' at the molecular bond 
length of 1.4 oq. As the density of deuterium increases to 
1.0 g/cm"^, one observes a drastic reduction in peak height 



TABLE I: FPEOS table with pressures and internal en- 
ergy per atom for deuterium. The statistical uncertain- 
ties from the PIMC simulations are in given in brack- 
ets e.g. 0.219(5)=0.219±0.005, 414.4(1. 6)=414.4±1. 6, or 
70230(400)=70230±400. 



Temperature 


Pressure 


Internal energy 


(K) 


(Mbar) 


(eV/atom) 



p = 1.96360xl0"^gcm-3 [r, = 14.0 ao] 

15625 0.001290(8) -10.83(3) 

31250 0.003364(8) -2.867(16) 

62500 0.009046(8) 11.930(13) 

95250 0.014640(8) 21.870(12) 

125000 0.019600(10) 30.080(15) 

181825 0.028920(8) 45.190(13) 

250000 0.040040(12) 63.04(2) 

400000 0.06440(2) 102.00(4) 

500000 0.08074(2) 128.10(4) 

1000000 0.16170(3) 257.30(6) 

2000000 0.32390(6) 515.90(10) 

4000000 0.64830(14) 1033.0(2) 

p = 3.11810xl0-^gcm-3 [r, = 12.0 ao] 

15625 0.002048(14) -10.97(3) 

31250 0.005105(14) -3.82(2) 

62500 0.013980(12) 10.860(12) 

95250 0.022970(11) 21.180(12) 

125000 0.030840(15) 29.490(15) 

181825 0.045690(13) 44.720(13) 

250000 0.063400(19) 62.67(2) 

400000 0.10220(4) 101.80(4) 

500000 0.12790(3) 127.70(4) 

1000000 0.25670(5) 257.10(5) 

2000000 0.51440(10) 516.00(10) 

4000000 1.0290(2) 1033.0(2) 

p = 5.38815xl0~^gcm-3 [r, = 10.0 ao] 

15625 0.00349(2) -11.280(20) 

31250 0.00845(2) -4.841(18) 

62500 0.02325(2) 9.401(14) 

95250 0.038870(19) 20.080(12) 

125000 0.05264(3) 28.620(16) 

181825 0.07841(2) 44.040(13) 

250000 0.10900(3) 62.090(19) 

400000 0.17630(7) 101.40(4) 

500000 0.22080(7) 127.30(4) 

1000000 0.44360(9) 257.00(5) 

2000000 0.88830(19) 515.50(11) 

4000000 1.7780(4) 1033.0(2) 
p = 1.05237xl0-^gcm-^ [r, = 8.0 ao] 

15625 0.00659(5) -11.550(17) 

31250 0.01580(6) -5.860(20) 

62500 0.04325(4) 7.477(13) 

95250 0.07371(4) 18.450(12) 

125000 0.10080(5) 27.220(16) 

181825 0.15150(4) 42.940(13) 

250000 0.21160(7) 61.18(2) 

400000 0.34300(13) 100.60(4) 

500000 0.43010(13) 126.70(4) 

1000000 0.86550(19) 256.50(6) 

?!Oonnnn i 7.3!sn('4~) ?,^?,M^(^^^ 
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TABLE II: TABLE I. {Continued.) TABLE III: TABLE I. (Continued.) 



Temperature 


Pressure 


Internal energy 


Temperature 


Pressure Internal energy 


(K) 


(Mbar) 


(eV/atom) 


(K) 


(Mbar) 


(eV/atom) 


p = 2.49451 xlO"^gcm- 


^ [r, = 6.0 ao] 


p = 0.199561 gcm"^ [r. 


= 3.0 Oo] 


15625 


0.01485(12) 


-11.820(17) 


15625 


0.124(3) 


-12.31(4) 


31250 


0.03556(11) 


-6.948(13) 


31250 


0.2740(17) 


-8.91(3) 


62500 


0.09550(10) 


4.773(12) 


62500 


0.6730(17) 


-1.11(3) 


95250 


0.16610(10) 


15.740(13) 


95250 


1.1760(16) 


8.14(3) 


125000 


0.23050(11) 


24.790(14) 


125000 


1.656(2) 


16.66(3) 


181825 


0.35170(12) 


40.920(16) 


181825 


2.6060(15) 


32.91(2) 


250000 


0.4946(2) 


59.40(3) 


250000 


3.753(2) 


51.98(3) 


400000 


0.8072(4) 


99.19(5) 


400000 


6.273(3) 


92.85(4) 


500000 


1.0150(3) 


125.40(4) 


500000 


7.943(2) 


119.60(4) 


1000000 


2.0480(5) 


255.50(6) 


1000000 


16.250(5) 


251.10(8) 


2000000 


4.1090(9) 


514.50(12) 


2000000 


32.760(7) 


510.90(11) 


4000000 


8.2310(17) 


1032.0(2) 


4000000 


65.740(14) 


1029.0(2) 


p = 4.31052 xlO"^gcm- 


^ [r, = 5.0 ao] 


p = 0.306563 gcm"^ [r. 


= 2.6 ao] 


15625 


0.02537(20) 


-11.970(16) 


15625 


0.219(5) 


-12.29(5) 


31250 


0.05980(19) 


-7.537(14) 


31250 


0.447(4) 


-9.15(4) 


62500 


0.15800(16) 


3.087(12) 


62500 


1.048(4) 


-1.90(4) 


95250 


0.27800(18) 


13.820(13) 


95250 


1.781(5) 


6.68(5) 


125000 


0.3880(3) 


22.910(18) 


125000 


2.509(4) 


14.98(4) 


181825 


0.5976(2) 


39.250(17) 


181825 


3.947(3) 


30.97(4) 


250000 


0.8460(3) 


58.01(3) 


250000 


5.693(5) 


49.91(5) 


400000 


1.3870(4) 


97.96(3) 


400000 


9.558(8) 


90.87(8) 


500000 


1.7460(5) 


124.40(4) 


500000 


12.120(6) 


117.70(6) 


1000000 


3.5310(8) 


254.50(6) 


1000000 


24.890(9) 


249.60(9) 


2000000 


7.0960(13) 


513.80(10) 


2000000 


50.250(13) 


509.70(13) 


4000000 


14.220(3) 


1031.0(2) 


4000000 


101.100(19) 


1030.0(2) 


p = 8.41898xl0"^gcni" 


^ [rs = 4.0 oo] 


p = 0.389768 gcm"^ [vs 


= 2.4 oo] 


15625 


0.0479(7) 


-12.21(2) 


15625 


0.298(12) 


-12.30(10) 


31250 


0.1150(5) 


-8.135(19) 


31250 


0.597(9) 


-9.21(7) 


62500 


0.2950(6) 


1.17(2) 


62500 


1.337(8) 


-2.40(7) 


95250 


0.5200(6) 


11.37(2) 


95250 


2.280(7) 


6.11(6) 


125000 


0.7308(7) 


20.29(3) 


125000 


3.175(8) 


14.09(7) 


181825 


1.1400(6) 


36.80(2) 


181825 


4.979(11) 


29.84(9) 


250000 


1.6250(9) 


55.71(3) 


250000 


7.206(9) 


48.84(7) 


400000 


2.6850(11) 


96.11(4) 


400000 


12.090(12) 


89.61(9) 


500000 


3.3910(11) 


122.70(4) 


500000 


15.370(14) 


116.70(11) 


1000000 


6.8810(15) 


253.30(6) 


1000000 


31.580(14) 


248.60(11) 


2000000 


13.850(3) 


513.00(11) 


2000000 


63.96(2) 


509.80(19) 


4000000 


27.750(6) 


1030.0(2) 


4000000 


128.40(3) 


1028.0(3) 


p = 0.1 


gcm"^ [rs = 


; 3.777 oo] 


p = 0.506024 gcirr-^ [r. 


= 2.2 Oo] 


15625 


0.0578(15) 


-12.21(5) 


15625 


0.42(3) 


-12.1(2) 


31250 


0.1351(6) 


-8.351(19) 


31250 


0.849(14) 


-9.16(9) 


62500 


0.3491(6) 


0.74(2) 


62500 


1.789(12) 


-2.68(7) 


95250 


0.6110(6) 


10.690(18) 


95250 


2.954(10) 


5.30(6) 


125000 


0.8598(8) 


19.57(2) 


125000 


4.088(12) 


13.02(7) 


181825 


1.3450(9) 


36.08(3) 


181825 


6.396(10) 


28.48(6) 


250000 


1.9200(9) 


55.00(3) 


250000 


9.243(12) 


47.20(7) 


400000 


3.1850(13) 


95.67(4) 


400000 


15.620(14) 


88.25(9) 


500000 


4.0180(11) 


122.10(3) 


500000 


19.84(3) 


115.10(18) 


1000000 


8.1680(17) 


252.90(5) 


1000000 


40.94(3) 


247.5(2) 


2000000 


16.440(3) 


512.40(11) 


2000000 


82.93(3) 


508.6(2) 


4000000 


32.970(7) 


1030.0(2) 


4000000 


166.50(4) 


1026.0(3) 



TABLE IV: TABLE L (Continued.) TABLE V: TABLE L {Continued.) 



Temperature 


Pressure 


Internal energy 


Temperature 


Pressure Internal energy 


(K) 


(Mbar) 


(eV/atom) 


(K) 


(Mbar) 


(eV/atom) 


p = 0.673518 gcm-^' 


[rs = 2.0 ao] 


p = 1.15688 gcm"^ [rs 


= 1.67 ao] 


15625 


0.59(4) 


-12.02(18) 


15625 


1.67(11) 


-10.9(3) 


31250 


1.28(2) 


-8.96(11) 


31250 


2.78(5) 


-8.35(14) 


62500 


2.461(10) 


-3.01(5) 


62500 


4.82(6) 


-2.87(16) 


95250 


3.930(7) 


4.43(3) 


95250 


7.05(6) 


3.49(17) 


125000 


5.413(6) 


11.91(3) 


125000 


9.65(6) 


10.93(17) 


181825 


8.446(8) 


27.08(4) 


181825 


14.45(6) 


24.77(17) 


250000 


12.200(7) 


45.63(3) 


250000 


20.81(5) 


42.94(14) 


400000 


20.660(14) 


86.64(6) 


400000 


35.12(5) 


83.36(13) 


500000 


26.270(15) 


113.50(7) 


500000 


44.82(3) 


110.50(9) 


1000000 


54.300(20) 


245.90(9) 


1000000 


92.99(5) 


243.30(13) 


2000000 


110.20(3) 


507.20(16) 


2000000 


189.00(11) 


505.0(3) 


4000000 


221.60(6) 


1026.0(3) 


4000000 


380.30(16) 


1024.0(4) 


p = 0.837338 [ 


r-s = 1.86 ao] 


p = 1.31547 [r^ 


= 1.60 ao] 


15625 


0.97(6) 


-11.3(2) 


31250 


3.46(9) 


-7.9(2) 


31250 


1.71(2) 


-8.85(8) 


62500 


5.65(10) 


-2.8(2) 


62500 


3.17(4) 


-3.17(14) 


95250 


8.31(8) 


3.75(19) 


95250 


5.04(4) 


4.30(14) 


125000 


10.97(6) 


10.44(15) 


125000 


6.80(3) 


11.43(11) 


181850 


16.66(6) 


24.76(14) 


181825 


10.52(2) 


26.29(9) 


250000 


23.73(5) 


42.50(12) 


250000 


15.15(3) 


44.67(10) 


400000 


39.92(7) 


82.72(17) 


400000 


25.62(4) 


85.52(16) 


500000 


50.66(7) 


109.20(17) 


500000 


32.67(4) 


112.70(17) 


1000000 


105.50(8) 


242.30(20) 


1000000 


67.44(6) 


244.9(2) 


2000000 


214.20(7) 


502.90(16) 


2000000 


136.90(11) 


506.2(4) 


4000000 


431.80(11) 


1022.0(3) 


4000000 


275.50(8) 


1026.0(3) 


p = 1.59649 gcm"^ [r. 


= 1.50 ao] 


p = 1.0 


gcm"^ [r, = 


= 1.753 ao] 


31250 


4.67(13) 


-7.5(3) 


15625 


1.33(7) 


-11.0(2) 


62500 


7.24(12) 


-2.7(2) 


31250 


2.22(4) 


-8.67(12) 


95250 


10.53(13) 


3.9(2) 


62500 


3.92(4) 


-3.23(14) 


125000 


13.68(11) 


10.4(2) 


95250 


6.06(4) 


3.88(13) 


181850 


20.19(7) 


23.87(13) 


125000 


8.16(4) 


10.90(11) 


250000 


28.78(8) 


41.57(15) 


181825 


12.57(3) 


25.60(10) 


400000 


48.31(7) 


81.53(14) 


250000 


18.07(3) 


43.85(10) 


500000 


61.33(10) 


107.90(19) 


400000 


30.36(3) 


84.00(9) 


1000000 


128.20(11) 


241.8(2) 


500000 


38.81(3) 


111.30(10) 


2000000 


260.40(15) 


503.1(3) 


1000000 


80.40(4) 


243.90(13) 


4000000 


524.30(13) 


1022.0(3) 


2000000 


163.20(6) 


505.00(17) 


p = 1.96361 gcm-^* [r. 


= 1.40 ao] 


4000000 


328.80(8) 


1024.0(3) 


31250 


6.4(2) 


-7.0(4) 


p = 1.00537 gcm"^ [r 


s = 1.750 oo] 


62500 


9.69(15) 


-2.1(2) 


15625 


1.35(9) 


-10.9(3) 


95250 


13.66(14) 


4.3(2) 


31250 


2.23(3) 


-8.69(10) 


125000 


17.11(10) 


10.05(16) 


62500 


4.03(4) 


-2.97(14) 


181850 


25.29(6) 


23.68(9) 


95250 


6.03(5) 


3.67(16) 


250000 


35.66(9) 


41.00(14) 


125000 


8.27(5) 


11.09(14) 


400000 


59.23(9) 


80.18(15) 


181825 


12.63(3) 


25.56(9) 


500000 


75.41(10) 


106.90(15) 


250000 


18.17(4) 


43.82(13) 


1000000 


156.80(9) 


239.50(15) 


400000 


30.54(5) 


84.05(16) 


2000000 


319.49(16) 


501.2(3) 


500000 


38.99(8) 


111.2(3) 


4000000 


644.54(18) 


1021.0(3) 


1000000 


80.80(7) 


243.7(2) 








2000000 


164.20(11) 


505.3(3) 








4000000 


330.50(12) 


1024.0(4) 









TABLE VI: TABLE L (Continued.) TABLE VIL TABLE I. (Continued.) 



Temperature 


Pressure Internal energy 


Temperature Pressure Internal energy 


(K) 


(Mbar) 


(eV/atom) 


(K) 


(Mbar) 


(eV/atom) 


p = 2.45250 gcm"^ [r„ 


= 1.30 ao] 


p = 10.0000 gcm"^ [r^ 


0.81373 ao] 


62500 


12.61(17) 


-2.2(2) 


125000 


142(3) 


19.0(8) 


95250 


18.18(16) 


4.9(2) 


181850 


182(2) 


31.8(7) 


125000 


22.09(15) 


10.05(19) 


250000 


225.7(1.7) 


44.5(5) 


181850 


32.63(10) 


24.00(13) 


400000 


334(4) 


80.5(1.3) 


250000 


45.03(12) 


40.55(15) 


500000 


414.4(1.6) 


106.2(5) 


400000 


74.54(10) 


79.73(13) 


1000000 


802(2) 


230.5(6) 


500000 


93.88(11) 


105.30(14) 


2000000 


1596.0(1.6) 


483.3(5) 


1000000 


195.60(17) 


238.1(2) 


4000000 


3276.0(1.6) 


1013.0(5) 


2000000 


398.20(18) 


499.4(2) 


8000000 


6592(3) 


2054.0(8) 


4000000 


804.2(3) 


1020.0(4) 


p = 10. 


5237 gcm-3 [r. 


= 0.80 ao] 


p = 3.11814 gcm"^ [rs 


= 1.20 ao] 


125000 


153(7) 


19.8(2.0) 


62500 


18.0(4) 


-1.4(4) 


181850 


197(4) 


33.1(1.0) 


95250 


24.6(3) 


5.4(4) 


250000 


242(3) 


45.4(9) 


125000 


30.1(5) 


11.0(5) 


400000 


353(4) 


80.4(1.2) 


181850 


43.0(3) 


24.4(4) 


500000 


434(4) 


105.3(1.1) 


250000 


57.8(4) 


39.9(4) 


1000000 


846(3) 


231.0(8) 


400000 


94.9(2) 


78.5(2) 


2000000 


1681(2) 


483.2(7) 


500000 


120.10(14) 


104.70(14) 


4000000 


3447(6) 


1011.0(1.8) 


1000000 


248.4(2) 


236.7(3) 


8000000 


6929(3) 


2051.0(9) 


2000000 


504.5(6) 


496.7(6) 


p = 15. 


7089 gcm"^ [rs 


= 0.70 ao] 


4000000 


1021.0(5) 


1018.0(5) 


181850 


346(9) 


40.0(1.7) 


p = 4.04819 gcm"^ [rs 


= 1.10 ao] 


250000 


419(9) 


55.0(1.8) 


62500 


26.2(1.1) 


-0.1(8) 


400000 


575(4) 


87.1(8) 


95250 


34.4(5) 


6.1(4) 


500000 


684(3) 


109.3(6) 


125000 


41.9(6) 


12.1(5) 


1000000 


1293(4) 


233.5(7) 


181850 


58.8(6) 


25.3(5) 


2000000 


2515(4) 


481.4(9) 


250000 


75.8(4) 


39.0(3) 


4000000 


5149(5) 


1011.0(1.1) 


400000 


125.1(3) 


78.5(2) 


8000000 


10390(5) 


2060.0(1.1) 


500000 


156.8(2) 


103.80(16) 


p = 24. 


9451 gcm"^ [rs 


= 0.60 ao] 


1000000 


321.0(3) 


234.0(3) 


400000 


1037(12) 


98.4(1.4) 


2000000 


651.5(7) 


492.9(5) 


500000 


1208(17) 


121(2) 


4000000 


1327.0(8) 


1018.0(6) 


1000000 


2133(11) 


239.2(1.4) 


p = 5.38815 gcm"^ [rs 


= 1.00 ao] 


2000000 


4025(9) 


480.9(1.1) 


95250 


51.9(1.6) 


8.3(9) 


4000000 


8195(16) 


1009.0(2.0) 


125000 


61.1(8) 


13.7(5) 


8000000 


16200(17) 


2020(2) 


181850 


81.8(1.2) 


25.9(7) 


16000000 


32950(13) 


4124.0(1.6) 


250000 


105.4(8) 


40.0(5) 


p = 43. 


1052 gcm-^* [r. 


= 0.50 ao] 


400000 


169.1(1.3) 


78.2(7) 


400000 


2212(30) 


123(2) 


500000 


212.2(1.3) 


104.1(8) 


500000 


2523(20) 


146.5(1.8) 


1000000 


429.4(1.1) 


233.5(6) 


1000000 


4002(30) 


256.4(2.0) 


2000000 


867.4(1.2) 


491.5(7) 


2000000 


7162(18) 


490.2(1.3) 


4000000 


1768.0(1.0) 


1018.0(6) 


4000000 


14180(17) 


1006.0(1.3) 


p = 7.39115 gcm"^ [r. 


= 0.90 ao] 


8000000 


28390(30) 


2044(2) 


95250 


80(2) 


10.3(9) 


16000000 


56880(20) 


4118.0(1.7) 


125000 


92(2) 


15.5(9) 


32000000 


114000(40) 


8273(3) 


181850 


120(2) 


27.4(1.0) 


64000000 


227900(90) 


16540(7) 


250000 


153.4(1.3) 


41.9(6) 


p = 84. 


1898 gcm"^ [rs 


= 0.40 ao] 


400000 


236.6(1.5) 


78.2(6) 


1000000 


9169(50) 


298.4(2.0) 


500000 


297.2(1.2) 


104.5(5) 


2000000 


14950(80) 


517(3) 


1000000 


590.2(1.4) 


231.8(6) 


4000000 


27960(70) 


1007(3) 


2000000 


1183.0(1.6) 


486.5(7) 


8000000 


54980(40) 


2019.0(1.4) 


4000000 


2422.0(1.6) 


1015.0(7) 


16000000 


110600(70) 


4093(3) 



32000000 222600(120) 8262(4) 
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TABLE VIII: TABLE I. {Continued.) 


Temperature 


Pressure 


Internal-energy 


(K) 


(Mbar) 


(eV/atom) 


p = 199.561 gcm"^ [r^ 


= 0.30 ao] 


2000000 


41350(800) 


597(12) 


4000000 


70230(400) 


1056(6) 


8000000 


129900(500) 


2000(8) 


16000000 


263400(300) 


4104(5) 


32000000 


527000(400) 


8247(7) 


64000000 


1049000(400) 


16450(7) 


p = 673.518 gcm"^ [r. 


= 0.20 ao] 


4000000 


299900(4000) 


1322(16) 


8000000 


504200(3000) 


2281(15) 


16000000 


897500(1300) 


4121(6) 


32000000 


1783000(1700) 


8251(8) 


64000000 


3569000(1600) 


16560(7) 


p = 1596.49 gcm"^ [r. 


= 0.15 ao] 


4000000 


1071000(20000) 


2002(40) 


8000000 


1555000(20000) 


2961(40) 


16000000 


2342000(17000) 


4517(30) 


32000000 


4565000(16000) 


8890(30) 


64000000 


8523000(8000) 


16670(17) 




4 

r (Bohr) 



FIG. 3: (Color online) The pair correlation functions g{r) 
derived from PIMC calculations: (a) the ion-ion correlation 
gddir); (b) the ion-electron correlation gde{r); (c) the electron- 
electron correlation geep{r) for parallel spins; (d) the electron- 
electron correlation geea{r) for anti-parallel spins, with differ- 
ent densities at 15 625 K. 



which demonstrates the pressure-induced dis sociat ion of 
D2 molecules, confirming earlier PIMC result^^^*^. This 
interpretation is also supported by the reduction of peak 
at r = in the gdeir) function in Fig.jsjb). Furthermore 
the positive correlation between pair of electrons with 
anti-parallel spin in Fig. [sf^d) is also disappearing with 
increasing density since they are no longer bound into 
molecules. Fig. [sjc) shows that there is always a strong 
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FIG. 4: (Color online) Pair correlation functions similar to 
Fig.jsjbut at a higher temperature of 2.5 x 10^ K and densities 
from 1.0 to 15.709 g/cm^. 
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FIG. 5: (Color online) Pair correlation functions similar to 
Fig. [3] but at a higher temperature of 2 x 10^ K and densities 
from 24.945 to 199.56 g/cm^ 



repulsion between electrons with parallel spins because 
of the Pauli exclusion principle but they approach each 
other more at higher densities. 

Figs. |4] and [5] show the pair correlation functions 
for different densities at much higher temperatures of 
2.5 X 10^ K and 2 x 10'' K. At these temperatures, D2 
molecules have completely dissociated as indicated by the 
absence of the peak in the ion-ion correlation function. 
The attractive forces between pair of ions have disap- 
peared and repulsion now dominates their interactions. 
At higher densities, particles are "packed" more tightly 
and approach each other significantly more so that the 
g{r) rise up more steeply and reach the values of 0.5 as 
much smaller distances. 

In Fig. |6j we compare the pair correlation functions for 
the fixed density of 10 g/cm"^ for temperatures ranging 
from 1.25 x 10^ K to 2 x lO*" K. It is interesting to note 
there is relatively little variation between the three curves 
below the Fermi temperature of Tp = 8.8 x 10^ K but 
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- T = 1 25,000 K 

- T= 250,000 K 

- T= 500,000 K 

- T = 2,000,000 K 
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FIG. 6: (Color online) Pair correlation functions similar to 
Fig. [3] at a fixed density of fO.O g/cm^ for temperatures rang- 
ing from 1.25 X 10^ K to 2 X 10^ K. 
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FIG. 7: (Color online) The comparisons of pressure (a) and 
internal energy (b) as a function of density from the the 
FPEOS, SESAME, KerleyOS tables. The error bars indicate 
the 1(7 statistical uncertainty in the PIMC simulations. Re- 
sults were normalized to non-interacting gas of classical ions 
and fermionic electrons. 



they differ significantly at 2 x 10^ K. This is a manifes- 
tation of Fermi degeneracy effects in which the electrons 
approach the ground state for temperatures well below 
the Fermi temperature. Then much of the temperature 
dependence of the pair correlation functions disappears. 
For example the pair correlation functions of electrons 
with anti-parallel spins are almost identical for the two 
lowest temperatures of 1.25 x 10^ K and 2.5 x 10^ K but 
they differ substantially from results at well above Tp. 
When the temperature raises above Tp, Pauli exclusion 
effects are reduced, the electrons start to occupy a vari- 
ety of states, which then has a positive feedback on the 
mobility of the ions. 



V. COMPARISONS OF THE FPEOS TABLE 
WITH SESAME AND KERLEY03 MODELS 

In this section, we compare the pressures and internal 
energies in our FPEOS table with predictions from the 
well-known semi-analytical SESAME and KerleyOS EOS 
tables. To illustrate how much the system deviates from 
an ideal plasma, we have normalized both pressure and 
energy to their corresponding values [Eid and Pid] of non- 
interacting gas of classical ions and fermionic electrons. 
This removes most of the temperature dependence and 
emphasizes the effects of the Coulomb interaction, which 
leads to a reduction in pressure and energy below the 
non-interacting values in all cases. 

In Figs. [7[-|9j we plot the pressure and the internal 
energy as a function of density for diff eren t tem peratures 



ranging from 31250 to 4 x 10 K. Figs. 10 - 13 show them 



as function of temperature for different densities varying 
between 0.1 and 84.19 g/cnx^. 

In Fig. [7) we compare FPEOS, SESAME, and KerleyOS 
results at a comparatively low temperature of 31250 K. 
This is difScult regime to describe by chemical mod- 
els because the plasma consists of neutral species like 



molecules and atoms as well as charged particles such as 
ions and free electrons. The interaction between neutral 
and charged species is very difficult analytically while 
it poses no major challenge to first-principles simula- 
tions. As is shown by Fig. [Tfa), the SESAME EOS 
predicts overall higher pressures at low density (p < 0.3 
g/cm^) but then all three models come to agree with each 
other at higher densities. The improved KerleyOS table 
still showed some discrepancy at very low densities, even 
though some improvements to the ionization equilibrium 
model have been mad^^. 

Fig. shows that the internal energies predicted 

by SESAME and KerleyOS are overall lower than the 
FPEOS values. The higher the density, the more dis- 
crepancy there is. Again, this manifests the difficulty of 
chemical models at such plasma conditions. 

One expects the pressure and internal energy to ap- 
proach the values of a non-interacting gas in the low- 
density and the high-density limit. At low density, parti- 
cles are so far away from each other that the interaction 
effects become negligible. At high density, Pauli exclu- 
sion effects dominate over all other interactions and all 
thermodynamic function can be obtained from the ideal 
Fermi gas. Just at an intermediate density range which 
still spans several orders of magnitude, the Coulomb in- 
teraction matters and significant deviations for the ideal 
behavior are observed. 

For a higher temperature of 2.5 x 10^ K, the pres- 
sure and energy are compared in Fig. |8] The low-density 
deuterium at this temperature becomes fully ionized and 
can therefore be described by the Debye-Hiickel plasma 
modeP^, which is based on the self-consistent solution of 
the Poisson equation for a system of screened charges. 
The pressure and energy per particle (counting electrons 
and ions) can be explicitly expressed as: 
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FIG. 8: (Color online) Same as Fig. [7] but for a diflterent 
temperature 2.5 x IC K. The Debye-Hiickel model is also 
shown for comparison. 



with the particle number density n, the Boltzmann con- 
stant k},, and the Debye length \d = \/ k^T / ^Trne^ . 

We have added the Debye-Hiickel results to Figs. [SpS] 
In Fig. [8] one finds that the simple Debye-Hiickel model 
perfectly agrees with our PIMC calculations in the lower 
densities up to 0.1 g/cm^, where the improved KerleyOS 
EOS also gives very similar pressures and energies. On 
the other hand, the SESAME EOS overestimates both 
pressure and energy even at such low densities. 

Fig. [sjja) exposes an artificial cusp in pressure in Ker- 
leyOS EOS at densities of 1.5 — 4 g/cm^ while the internal 
energy curve is smooth. This artificial pressure cusp ap- 
pears for all temperatures at roughly the same density 
and may be related to the artificial double compression 
peaks in the principal Hugoniot predicted by KerleyOS 
E0&i2.. The Debye-Hiickel model fails at densities higher 
than 0.2 g/cm'^ for this temperature. It is only applica- 
ble to weakly interacting plasmas but otherwise predicts 
unphysically low pressures and energies. 

As the temperature increased to 4 x 10^ K, the Debye- 
Hiickel model agrees very well with FPEOS in both pres- 
sure and energy over a wide range of densities up to 20 
g/cm'^ as shown in Fig. [9] Significant differences in both 
pressure and energy are again found for the SESAME 
EOS, when compared to FPEOS and KerleyOS tables. It 
should also be noted that the internal energy predicted 
by KerleyOS is slightly lower than those of FPEOS and 
the Debye-Hiickel model for p = 0.1 — 20 g/cm'^. 

In Figs. [T0]-[T3] we compare the pressure and energy 
versus temperature for specific densities of 0.1, 1.0, 10.0, 
84.19 g/cm'^. At high temperature where the plasma is 
fully ionized, the Debye-Hiickel model well reproduces 
the FPEOS pressures and energies very well. It is in- 
teresting to note that the SESAME table overestimates 
the pressure and energy even for a fully ionized pl asma at 
densities greater than 1.0 g/cm"^ as shown in Figs 
For a very low density of 0.1 g/cm"^ 



Fig. 
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shows that 



the improved KerleyOS agrees very well with FPEOS, 
while the SESAME results are noticeably higher. More- 
over, the improvements made to KerleyOS have resulted 
in remarkable agreement with FPEOS for intermediate 
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FIG. 9: (Color online) Same as Fig. [7] but for a different 
temperature of 4 x 10® K. 
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FIG. 10: (Color online) Pressure (a) and energy (b) as a func- 
tion of temperature from FPEOS, SESAME, and KerleyOS 
tables for deuterium density of 0.1 g/cm^. 



densities of 0.1 and 10.0 g/cm^ depicted by Figs. 11 and 
|12| Only a small deviation in the internal energy be- 
tween KerleyOS and our FPEOS results can be found at 
the lowest temperature for 1.0 g/cm^. 

At a higher density of 84.19 g/cm^, the SESAME EOS 
again significantly deviates from both the FPEOS and 
the KerleyOS EOS as is illustrated by Fig. [13| The lat- 
ter two EOS tables give very similar results in internal 
energy almost for the entire temperature range, though 
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FIG. 11: (Color online) Same as Fig. 10 but for a different 
deuterium density of 1.0 g/cm^. 
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FIG. 12: (Color online) Same as Fig. 10 but for a higher 
deuterium density of 10.0 g/cm^. 
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FIG. 13: (Color online) Same as Fig. 10 but for much higher 
deuterium density 84.19 g/cm^. 



the pressures predicted by KerleyOS are higher than the 
FPEOS ones for temperatures varying from 2 x 10^ K to 
2 X 10^ K. In contrast to the significant EOS differences 
seen from SESAME, the improved KerleyOS table is over- 
all in better agreement with the FPEOS table, although 
subtle discrepancies and an artificial pressure cusp still 
exist in the KerleyOS EOS. 



VI. APPLICATIONS TO ICF 

With the EOS comparisons discussed above, we now 
investigate what differences can be observed when these 
EOS tables are applied to simulate ICF shock timing 
experiments and target implosions. Using radiation- 
hydrodynamics codes, both one-dimensional LILAC^ 
and two-dimensional DRACC^, for simulations of ex- 
periments, we can explore the implications of our first- 
principles equation of state table for the understanding 
and design of ICF targets. 

We first study the shock timing experiment d^^ l ^^ * per- 
formed on the Omega laser facility. As the fuel entropy in 
ICF implosions is set by a sequence of shocks, the timing 
of shock waves in liquid deuterium is extremely impor- 
tant for the ICF target performance. In shock timing 
experiments, the carbon deuterium (CD) spherical shell. 



FIG. 14: (Color online) (a) The triple-picket pulse shape for 
shock-timing experiments using cryogenic deuterium; and (b) 
the measured shock speed in liquid-D2 (black solid-line) com- 
paring with hydro-dynamics predictions using the SESAME, 
FPEOS, KerleyOS, and QEOS models. 



900 /im in diameter and 10 /im thick, in a cone-in-shell 
geometrj*^ were filled with liquid deuterium. VISAR 
(velocity interferometery system for any refiector) was 
used to measure the shock velocity. As is shown in 
Fig. [I4j[a) , the triple-picket laser pulses are designed to 
launch three shocks into the liquid deuterium. The ex- 
perimental results are plotted in Fig. [l4|b) , in which the 
shock front velocity is shown as a function of time. One 
finds that when the second shock catches up the first 
one at around 1.5 ns, the shock-front velocity exhibits 
a sudden jump. Another velocity jump at 2.2 ns oc- 
curs when the third strong shock overtakes the previous 
two. With the hydro-code LILAC, we have simulated the 
shock timing experiments using different EOS tables in- 
cluding FPEOS, SESAME, KerleyOS, and QEO^. The 
radiation hydro-simulations have used the standard flux- 
limited (/ = 0.06) thermal transport model, although 
a nonlocal model has resulted in better agreement with 
experiment for the speed of first shoclP^. The results of 
FPEOS, SESAME and KerleyOS are in good agreement 
with the experimental observation, while the QEOS pre- 
dicts s much lower shock velocity and early catching up 
time. The shock timing experiments can only explore a 
small range of deuterium densities (0.6—2.5 g/cm"^) and 
temperatures (3 — 10 eV). In these plasma conditions the 
SESAME and KerleyOS have been adjustecPl to match 
to the first-principles calculations, which can be seen in 
Fig. [Tl] Thus, the shock velocity differences predicted 
by the FPEOS, SESAME, and KerleyOS are very small 
in such plasma conditions. 

Next, we examine the implications of coupling and 
degeneracy effects in ICF implosions. The possi- 
ble differences in target compression and fusion yields 
of ICF implosions are investigated through radiation 
hydro-simulations using FPEOS in comparison to re- 
sults predicted by SESAME and KerleyOS. The LILAC- 
simulation results are compared in Figs. [15] and |16[ re- 
spectively, for a DT implosion on OMEGA and a hydro- 
equivalent direct-drive design on the NIF. In Figs.fl5Fa) 



12 



and [l6|[a), we plot the laser pulse shapes consisting of 
triple-pickets and the step-main-pulse. The cryogenic 
OMEGA DT target (860 fim diameter) has a 10-^m 
deuterated plastic ablator and ^ 65 fim of DT ice. Fig- 
ure 15 b) shows the density and temperature profiles at 
the end of the laser pulse {t — 3.8 ns) from the FPEOS 
(red/solid line), the SESAME (green/dashed line), and 
the KerleyOS (blue/dot-dashed line) simulations. AT this 
time the shell has converged to a radius of ~ 160 from 
its original radius of ~ 430 fim. The shell?s peak den- 
sity and average temperature were p ~ 5.6 g/cm^ and 
T ~ 21 eV, which correspond to the coupled and degen- 
erate regimes with F ~ 1.22 and 9 ~ 0.47. It is shown 
that the FPEOS simulation predicted ~ 10% lower peak 
density but 15% higher temperature relative to the 
SESAME prediction. As is shown by the comparisons 
made in Fig. [s] and in Ref. pj, the FPEOS predicts 
slightly stiffer deuterium than SESAME at the similar 
temperature regime. This explains the lower peak den- 
sity seen in Fig. 15 'b). The ~ 15% higher tempera- 
ture in the FPEOS case was originated from the lower 
internal energy [see Fig. |8]^b)]. Since the laser abla- 
tion does the same work/energy to the shell compres- 
sion and its kinetic motion, a lower internal energy in 
FPEOS means more energy is partitioned to heat the 
shell, thereby resulting in a higher temperature. Such a 
temperature increase and density drop can have conse- 
quences in the implosion performance. Despite the sub- 
tle EOS differences discussed above, the KerleyOS simula- 
tion show very similar results when compared to FPEOS. 
Only small differences in temperature profile can be seen 
between the FPEOS and KerleyOS simulations, both of 
which are in remarkable contrast to the SESAME case. 
Figure 15 c) show the density profile at the peak com- 



pression, in which the predicted peak density (pp ~ 210 
g/cm'^) is ^ 25% lower according to FPEOS and Ker- 
leyOS compared to the SESAME prediction {pp ~ 260 
g/cm'^). The history of areal density pi?-evolution and 
neutron production were shown in Fig. 15 'd). One sees 



that the peak pR and neutron yield are also reduced by 
- 10% - 20% when the FPEOS and KerleyOS are com- 
pared to the SESAME predictions. The absolute neutron 
yield drops from 8.44 x 10" predicted by SESAME to 
~ 6.91 X 1013 (FPEOS) and - 6.93 x lO^^ {KerleyOS). 

Figure [16] shows the similar effects for the hydro- 
equivalent direct-drive NIF design with 1-MJ laser en- 
ergy. The NIF target (cj) — 2.954-mm) consists of 27-pm 
plastic ablator and 170-pm DT ice. The triple-picket 
drive pulse has a total duration of '-^ 11.4 ns and a peak 
power of ~ 240- TW. We also found a decrease in pp and 
a slight temperature increase for the FPEOS and Ker- 
leyOS relative to SESAME simulations near the end of the 
laser pulse (t — 9.2 ns), shown by Fig. [T6j[b). The peak 
density at the stagnation dropped from 481 (SESAME) 
to ~ 445 g/cm^ (FPEOS/Keriey03), which is indicated 
by Fig. 16 c). The resulting pR and neutron yield as 

The yield 
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FIG. 15: (Color online) The hydro-code simulations of a cryo- 
genic DT implosion on OMEGA using the three different EOS 
tables including SESAME, FPEOS, and KerleyOS: (a) The 
laser pulse shape; (b) The density-temperature profiles of the 
imploding DT shell at the middle of main laser pulse (i = 3.8 
ns); (c) the density profile at the peak compression (t=4.26 
ns); and (d) the areal density pR and yield as a function of 
time. 
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FIG. 16: (Color online) Similar to Fig. [Ts] but for a hydro- 
equivalent direct drive, 1 MJ ignition design for NIF. 



Y ~ 1.57 X 10^9 (FPEOS) and Y ~ 1.55 x lO^^ (Ker- 
leyOS). Consequently, the energy gain decreased from 
49.1 (SESAME) to 44.2 (FPEOS) and 43.8 (KerleyOS). 
It is noted that the ^11% gain reduction for this design 
is much modest than the 1.5-MJ NIF design discussed in 
Ref. [^"] in which more than ^ 20% gain difference has 
been seen between FPEOS and SESAME simulations. 
This is attributed to the different density-temperature 
trajectories that the two designed implosions undergo, in 
which the EOS variations among FPEOS, SESAME and 
KerleyOS are different. 
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FIG. 17: (Color online) The density contour plots at t=3.85 
ns from two-dimensional DRACO simulations of the OMEGA 
cryogenic DT-implosion shown in Fig. |15[ respectively using 
the SESAME (a) and the FPEOS (b) for the fuel DT. It 
is noted that the various perturbation sources have been in- 
cluded up to mode I = 150. 



E 




Mode f. Time (ns) 

FIG. 18: (Color online) (a) The modal spectrum of pertur- 
bation amplitude at the ablation surface for t ~ 3.0 ns (the 
beginning of acceleration); (b) The perturbation growth in 
pR as a function of time, analyzed from the 2D-DRACO sim- 
ulations shown in Fig. [T7| 



Finally, we discuss the implications of the coupling and 
degeneracy effects in FPEOS to IGF target performance 
beyond the ID physics studied above. As we knew that 
various perturbations seeded by target roughness and 
lasers can grow via the Rayleigh-Taylor (RT) instabilitjJ^ 
during the shell acceleration/deceleration phases in ICF 
implosions, it is important to properly simulate the RT 
growth of fusion fuel for understandin g tar get perfor- 
mance (compression and neutron yieldsp^'^. Since the 
RT growth depends on the compressibility of materi- 
als, the accurate equation-of-state of deuterium is essen- 
tial to ICF designs. As an example, we have used our 
two-dimensional radiative hydro-code DRACO to simu- 
late the cryogenic DT implosion on OMEGA [discussed 
in Fig. 15 . The various perturbation sources, includ- 



ing the target offset, ice roughness, and laser irradiation 
non-uniformities measured from experiments, have been 
taken into account up to a maximum mode of I = 150. 
We have compared the FPEOS and SESAME simulation 
results in Fig. 17 for t = 3.85 ns near the end of ac- 



celeration, in which the density contours are plotted in 
the Zi?-plane (azimuthal symmetry with respect to the 
Z-axis is assumed). Visible differences in the DT shell 
density can even be seen by eye from Figs.[T7|^a) and (b). 
The FPEOS simulation resulted in more "holes" and den- 
sity modulations along the shell than the SESAME case. 

To further analyze the perturbation amplitudes, we 
have decomposed the ablation-surface modulations into a 
modal spectrum that is shown in Fig. 18 'a), at the start 
of sheU acceleration (t=3.0 ns). We find that the FPEOS 
predicted larger amplitudes than the SESAME case al- 
most over the entire modal range. As the deuterium 
Hugoniot was shown in Ref. p^, the FPEOS predicted 
softer deuterium compared to SESAME for pressures be- 
low ^ 2 Mbar. Thus, the softer deuterium can be more 
easily "imprinted" by the perturbations brought in via 
the series of shocks. This results in larger perturbation 
amplitudes in FPEOS than SESAME simulations. The 
Rayleigh-Taylor instability further amplifies these per- 



turbations during the shell acceleration. As indicated by 
Fig. 18 ^b), the arms of fuel pR modulation increases to 
a few mg/cm^ at the end of the laser pulse. These per- 
turbations penetrated into the inner surface of the DT- 
shell will become the seeds for further RT growth during 
the shell's deceleration phase. They eventually distort 
the hot-spot temperature and density, thereby reducing 
the neutron production. At the end, we found that the 
SESAME simulation resulted in a neutron- averaged ion 
temperature of (T^) — 2.6 keV and a neutron yield of 
y = 5.2 X 10^^; while due to the larger perturbations pre- 
dicted the FPEOS simulation has given an (T^) = 2.3 keV 
and neutron yield oiY = 3.7 x 10^^, which is more close 
to experimental observations of (Tj) = 1.8 ±0.5 keV and 
y = 1.9 X 10^^ 



VII. SUMMARY 

In conclusion, we have derived a first-principles equa- 
tion of state table of deuterium for ICF applications from 
PIMC calculations. The derived FPEOS table covers 
the whole plasma density and temperature conditions 
in low-adiabat IGF implosions. In comparison with the 
chemical model based SESAME table, the FPEOS table 
show significant difference in internal energy and pres- 
sure for coupled and degenerate plasma conditions; while 
the recently improved KerleyOS table exhibited fewer and 
smaller discrepancies when compared to the FPEOS pre- 
dictions temperature higher than ^ 10-eV. Although 
subtle differences at lower temperatures (T < 10 eV) and 
moderate densities (1 — 10 g/cm'^) have been identified 
and an artificial pressure cusp still exists in the KerleyOS 
table, radiation hydro-simulations of cryogenic IGF im- 
plosions using the FPEOS and KerleyOS tables have given 
similar peak density, areal density pR, and neutron yield, 
which are remarkably different from the SESAME sim- 
ulations. Both the FPEOS and the KerleyOS predicted 
~ 25% less peak density, ~ 10% smaller pR, and ~ 10%- 
20% less neutron yield, when compared to the SESAME 
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case. Two-dimensional simulations further demonstrated 
the significant differences in target performance between 
the FPEOS and SESAME simulations. In general, the 
FPEOS simulations resulted in better agreement with 
experimental observations in terms of ion temperature 
and neutron yield. It is also noted that the extreme 
conditions covered by the FPEOS table are also impor- 
tant in astrophysics and planetary sciences, for example, 
to model the evolution of staraSHl and to understand the 
thermodynamical properties of stellar mattei'^. 
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